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Abstract. Dust depletion plays a key role in understand- 
ing the nature of Damped Lyman a systems(DLAs). In 
this paper we point out a previously unnoticed anticor- 
relation between the observed abundance ratio [X/Zn] 
(where Zn is assumed to be undepleted and X stands for 
the refractories Fe, Cr and Ni) and metal column density 
([Zn/H]+log(Nff/)) in DLAs. We suggest that this trend is 
an unambiguous sign of dust depletion, since metal column 
density is a measure of the amount of dust along the line of 
sight. Assuming that DLAs are (proto-) galactic disks and 
using detailed chemical evolution models with metallicity 
dependent yields we study chemical evolution and dust de- 
pletion patterns for a and iron— peak elements in DLAs. 
When observational constraints on the metal column den- 
sity of DLAs are taken into account (as suggested in Boisse 
et al. 1998) we find that our models reproduce fairly well 
the observed mild rcdshift evolution of the abundances of 
8 elements (Al, Si, S, Cr, Mn, Fe, Zn and Ni) as well as 
the observed scatter at a given redshift. By considering the 
aforementioned dependence of abundance ratios on metal 
column density, we further explore the general dust deple- 
tion pattern in DLAs, comparing to our model results and 
to a solar reference pattern. We find that for low metal col- 
umn densities (no depletion), our models compare fairly 
well to the data, while a solar pattern has difficulties with 
Mn. At high metal column densities (amount of deple- 
tion ~0.5 dex), the solar pattern describes the data quite 
well, while our models have difficulties with S. We suggest 
that further measurements of those key elements, i.e. Zn, 
S and Mn, will help to gain more insight into the nature of 
DLAs. In any case, the presently uncertain nucleosynthe- 
sis of Zn in massive stars (on which a large part of these 
conclusions is based) should be carefully scrutinised. 

Key words: ISM: abundances — ISM: dust — Galaxies: 
spirals — Galaxies: abundances — Galaxies: evolution 



1. Introduction 

Damped Lyman a systems (DLAs) are high column 
density gaseous systems (Nhi>2 10 20 cm -2 ), detected 
through their absorption lines in the optical spectra of 
quasars, up to relatively high redshifts (up to z *~5). 
Their study constitutes a powerful means to investigate 
the properties of distant galaxies (or of their building 
blocks). In particular, DLA metal abundances have been 
widely used in the past few years to probe the nature of 
DLAs (e.g. Lauroesch et al. 1996; Lu et al. 1996; Pettini ct 
al. 1994, 1997b, 1999, 2000; Prochaska & Wolfe 1999,2000, 
and many others)). However, it is not clear whether the 
observed abundances allow that, because of various biases: 
depletion of metals into dust (Pei and Fall 1995) or bias 
against too high or too low metal column densities (Boisse 
et al. 1998). The latter bias, in particular, may explain the 
observed absence of evolution in the absolute abundance 
of Zn/H as a function of redshift, as suggested in Prantzos 
and Boissier (2000a, hereafter PB2000a). 

Abundance ratios offer a better diagnostic tool than 
absolute abundances for the study of the chemical evo- 
lution of a system. However, in the case of DLAs the 
possibility of depletion into dust grains (even to small ex- 
tent) makes difficult a direct interpretation of the observed 
abundance patterns. Several attempts have been made in 
recent years to recover the intrinsic abundance pattern 
from the observed one and to compare the result to well 
known patterns of stars in the Milky Way (e.g. Vladilo 
1998, 1999; Pettini et al. 2000), but with rather modest 
success up to now. 

In this work we reassess the question of abundance 
patterns in DLAs. First, we point out an anticorrelation 
between the observed abundance ratio X/Zn (where Zn is 
assumed to be undepleted and X stands for the refracto- 
ries Fe, Cr and Ni) and metal column density (Sec. 2). We 
suggest that this, previously unnoticed, trend is an unam- 
biguous sign of depletion, since metal column density is a 
measure of the amount of dust along the line of sight. We 
then study the chemical evolution of those systems, as- 
suming they are (proto-) galactic disks. This assumption 
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is rather controversial, since other systems (e.g. Low Sur- 
face Brightness galaxies, dwarf irregulars, galactic haloes; 
see Jimenez et al. 1998, Matteucci et al. 1997, Valageas et 
al. 1999, Haehnelt et al. 1998, Petitjean and Ledoux 1999) 
have also been suggested; however, the various zones of our 
multi-zone disk models have sufficiently diverse chemical 
histories to account for systems other than the Milky Way, 
a "prototype" disk galaxy. Our models, briefly presented 
in Sec. 3.1, have been shown to reproduce successfully a 
large number of observed properties of spirals at high and 
low redshifts. An important ingredient of our study is the 
use of metallicity dependent yields for massive stars and 
SNIa (presented in Sec. 3.2). In Sec. 4 we show that our 
results reproduce fairly well both the observed slow evo- 
lution with redshift of the absolute abundances of several 
elements, and the corresponding scatter; observational bi- 
ases, as already suggested in PB2000a for the case of Zn, is 
a key to that success. In Sec. 5 we compare our results to 
observed abundance ratios in DLAs, taking into account 
the aforementioned trend of depletion vs. metal column 
density. We find that for low metal column densities (i.e. 
no depletion) our models compare fairly well to the data, 
while a solar pattern has difficulties with Mn. At high col- 
umn densities (small depletion) the solar pattern describes 
the data quite well, while our models have difficulties with 
S. We suggest that further measurements of those key el- 
ements, i.e. Zn, S and Mn, will help to gain some insight 
into the nature of DLAs. In any case, the presently uncer- 
tain nucleosynthesis of Zn (on which a large part of these 
conclusions is based) should be carefully scrutinised. 

2. DLA abundance ratios vs. metal column 
density 

Studies of the elemental abundances of the interstellar 
medium (ISM) in the Milky Way have shown that the 
gas phase abundances of most elements are below the cor- 
responding solar values. This underabundance is generally 
attributed to depletion into dust grains. The amount of de- 
pletion depends on the nature of the interstellar medium. 
In cold, dense, disk gas, Mn, Cr, Fe and Ni are consid- 
erably depleted (by 1.5 to 2 dcx), while S/H is nearly 
solar (Savage and Sembach 1996 and references therein). 
In warm disk gas, the amount of depletion is smaller (~ 
1 dex for Fe— peak elements), while in warm, halo gas, a 
negligible amount of depletion is found (<0.5 dex) for all 
elements, except Fe and Ni which are depleted by 0.6 to 
0.8 dex. 

A more quantitative assessement of the amount of de- 
pletion can be made by plotting the abundance of each 
ionic species as a function of the hydrogen column density 
Nhi along the line of sight. In a recent work, Wakker and 
Mathis (2000) compiled data from several studies concern- 
ing high column densities (Nhi = 10 20 to 7 10 21 cm~ 2 ) 
and extended that range down to 10 18 cm~ 2 by using 
data from high— and intermediate— velocity gas clouds 



(most of them have intrinsically near-solar abundance). 
They showed that over the whole range of HI column den- 
sities (covering almost three orders of magnitude), there 
is a tight anticorrclation between log(X/H) and log(Nni) 
for several elements: high column densities correspond to 
a larger amount of depiction, with a remarkably low dis- 
persion around each value of Nhi • The slope of the anti- 
correlation depends on the ionic species considered; it is 
-0.39±0.04 for Mnll and -0.59±0.04 for Fell, two species 
that have also been measured in DLAs. 

It is interesting to see whether such an anticorrelation 
exists for DLAs. In Table 1 we present a compilation of 
currently available data for DLA abundances that we use 
in this work. The redshift range of the DLAs and their 
number in each survey are also given in the Table. Data 
displayed in all our figures are directly taken from the 
source references, with no attempt to homogenize them. 



Table 1. a and Fe— peak element abundances in DLAs 



Ref. 


Zabs 


obj. 


Al 


Si 


s 


Cr 


Mn 


Fc 


Ni 


Zn 


(1) 


0.3 - 1.0 


6 
















X 


(2) 


1.7 - 4.4 


15 




X 


X 






X 




X 


(3) 


1.9 - 2.5 


6 




X 


X 


X 




X 




X 


(4) 


0.68, 1.15 


2 




X 






X 


X 




X 


(5) 


0.7 - 4.5 


23 


X 


X 


X 


X 


X 


X 


X 


X 


(6) 


2.3090 


1 






X 






X 






(7) 


3.3901 


1 




X 




X 




X 


X 


X 


(8) 


1.7 - 3.0 


15 








X 








X 


(9) 


0.7 - 3.3 


17 








X 








X 


(10) 


1.1 - 1.5 


5 




X 




X 


X 


X 


X 


X 


(11) 


0.6 - 1.5 


6 




X 




X 


X 


X 


X 


X 


(12) 


1.7 - 4.2 


22 


X 


X 


X 


X 




X 


X 


X 


(13) 


2.4 - 4.2 


4 












X 







References: (1) Boisse et al. 1998; (2) Centurion et al. 1998; 
(3) Centurion ct al. 2000; (4) dc la Varga et al. 2000; (5) Lu 
et al. 1996; (6) Molaro et al. 1998; (7) Molaro et al. 2000; (8) 
Pettini et al. 1994; (9) Pettini et al. 1997b; (10) Pettini et al. 
1999; (11) Pettini et al. 2000; (12) Prochaska & Wolfe 1999; 
(13) Prochaska & Wolfe 2000 

In order to evaluate the depletion in DLAs as a func- 
tion of column density, one cannot use the absolute abun- 
dances [X/H] because DLAs have different metallicities: a 
low [X/H] value may correspond either to large depletion 
or to little chemical evolution (or both) . For those reasons, 
we adopt here the [X/Zn] ratio as a measure of depletion 
of a given element X, assuming that Zn is always unde- 
plctcd, irrespectively of the HI column density. By using 
[X/Zn] to measure depletion we also assume that the real 
[X/Zn] value of each DLA (i.e. before depletion) is always 
the same, which is probably not the case. However, we 
are looking here for systematic effects as a function of col- 
umn density; any deviations of the real [X/Zn] from the 
value adopted here ([X/Zn] = 0) would merely manifest 
themselves as dispersion around an average trend. 

On the left panels of Fig. 1 we present the DLA data 
of [X/Zn] as a function of Nhi column density. No anticor- 
relation is found, contrary to the situation in gas clouds 
in our Galaxy, which is shown by solid lines in the panels 



Please give a shorter version with: \authorrunning and/or \titilerunning prior to \maketitle 



3 



of Mn and Fe (from the fits of Wakker and Mathis 2000). 
The absence of such an anticorrclation in DLAs, compared 
to local measurements, seems puzzling. Does it mean that 
depiction (if any) in DLAs is independent of HI column 
density? 

Here, we suggest that a proper measurement of de- 
pletion is provided not by the HI column density, but 
rather by the column density of (undepleted) metals, e.g. 
by log(N Zn ) = [Zn/H]-7.83+log(N ffl ), where -7.83 = 
log(Zn/H) Q . In the case of clouds in our Galaxy, metal- 
licity is always the same ([Zn/H] ~ 0), making it possible 
to use Nhi instead of Nz n - However, in DLAs, because of 
the different [Zn/H], Nhi is not proportional to Nz n and 
cannot be used to measure depletion. 

When [X/Zn] is plotted as a function of T = 
[Zn/H]+log(Nni) (right panels in Fig. 1), one can see that: 
i) there is no dependence of Si/Zn and S/Zn on T ', as ex- 
pected for those two elements which are also found not to 
be substantially depleted in local gas clouds, and ii) there 
is a clear anticorrelation between [X/Zn] and T for the re- 
fractory elements Cr, Fe and Ni (as can be seen from the 
dashed lines) while the scarce data for Mn do not allow 
for a conclusion. 

We think that the anticorrelation obtained between 
[X/Zn] and T in DLAs is an unmistakable sign of depletion 
for Fe, Cr and Ni. Metal column density is the appropriate 
factor, since it measures the total amount of metals along 
the line of sight, and thus the amount of dust; Nhi alone 
is inappropriate, since further assumptions about the dust 
to gas ratio have to be made. 

How can one explain the fact that [X/Zn] seems to be 
independent of Nhi for all elements in DLAs? We note that 
there exists a well known anticorrelation between metal- 
licity, as expressed by [Zn/H], and HI column density in 
DLAs. This anticorrelation, first noticed by Boisse et al. 
(1998), is attributed to selection biases (since Zn is not de- 
pleted into dust grains): regions of high [Zn/H] and high 
Nhi are opaque to light from background quasars and 
thus undetectable, while regions of low [Zn/H] and low 
Nhi (i.e. of low Zn column density) are not detectable by 
current DLA surveys. In PB2000a we substantiated this 
argument by running detailed chemical evolution models 
for disk galaxies. We showed how this "filter" can lead to 
a "no-evolution" picture for the average DLA metallicity 
as a function of redshift, in agreement with observations. 
We shall return to that question in Sec. 3.3, presenting 
results for other elements than Zn. Here we simply note 
that, since high Nhi DLAs are asociated, on average, with 
low [Zn/H] and vice— versa, the depleted ratio [X/Zn] is 
expected to be, on average, independent of Nhi, as ob- 
served. On the contrary, when [X/Zn] is plotted as a func- 
tion of T (which is directly proportional to the amount of 
dust), the anticorrelation expected by local observations 
becomes clear. 

An inspection of the Fe panel in Fig. 1 shows that 
the absolute [Fe/Zn] values are systematically higher in 



DLAs than in the Milky Way clouds (by ~ 0.7 dex); a 
similar remark also holds for Mn, although the scarcity 
of data makes its case less compelling. If metal column 
density is the relevant factor, as argued in the previous 
paragraphs, how can this systematic offset (by factor of <~ 
5) be understood? We see two possibilities: 

1) The undepleted [Fc/Zn] in DLAs is not zero (i.e. the 
Fc/Zn ratio is not solar) as assumed here, but ~ 0.7 (i.e. 5 
times solar); however, despite current uncertainties on the 
nucleosynthesis of Zn, this possibility is quite implausible. 
Indeed, Fe/Zn is solar on average in both halo and disk 
stars of the Milky Way. 

2) Wakker and Mathis (2000) suggest that the very 
small scatter around the average depletion pattern in the 
Milky Way clouds (within dotted lines in the Fe panels of 
Fig. 1) is due to the fact that one single cloud dominates 
the column density of each line of sight. In that case, a 
single value of (depleted) Fe/Zn corresponds to a given 
value of Nhi (or of T, since [Zn/H] ~ solar in the Milky 
Way clouds). In the case of DLAs, it is tempting to assume 
that a given line of sight intercepts several clouds: the 
depletion of [Fe/Zn] within each cloud would be small, but 
the total column density of metals (!F) would be high. This 
picture could explain the offset of DLAs with respect to 
the local pattern in Fig. 1. Taking into account the error 
bars in [Fe/Zn] values and the possibility of systematic 
errors, it is difficult to tell whether a couple of clouds or 
a larger number is necessary. 

At this point, we notice that in the case that DLAs 
are composed by many systems of individual column den- 
sities logN(HI)<19.5, ionisation effects could play an im- 
portant role in shaping the observed abundance ratios. 
As suggested recently by Howk and Sembach (1999), the 
observed ratios may systematically overestimate the in- 
trinsic ones in the cases of Si/Fe and Zn/Cr, if ionised gas 
is present in these systems, "thereby mimicking the effects 
of alpha-element enrichment or dust depletion" . Izotov et 
al. (2000) have further explored this issue, on the basis 
of similarities observed between some DLAs and nearby 
blue compact dwarf galaxies; using a schematic two zone 
model, they found that ionisation effects may play some 
role in shaping the observed abundances of these systems, 
although they recognised that their model suffered from 
various uncertainties related to the unknown structure and 
geometry of the absorbing systems, the spatial distribu- 
tion of the ionising sources and their radiation field. We 
do not consider any further this possibility here, since we 
think that the anticorrelation displayed in Fig. 1, i.e. the 
dependence of X/Zn on Zn column density for well known 
refractory elements (like Fe, Ni and Cr) but not for e.g. 
Si, suggests that dust depletion is the main factor. 

The observations of Fig. 1 tell nothing about the real 
(undepleted) abundances in DLAs; further assumptions 
have to be made for that. The investigation of the abun- 
dances and the depletion patterns in DLAs by means of a 
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Fig. 1. Gas phase abundance ratios [X/Zn] in DLAs, plotted as a function of HI column density (left panels) and of 
Zn column density (right panels). Data are from references given in Table 1 and corresponding typical uncertainties 
are shown in each one of the right panels by error bars. No correlation is seen in the left panels, while in the right 
panels an anticorrelation is found for the refractory elements Cr, Fe and Ni (shown by dashed lines, obtained by a least 
squares fit to the data) . In the panels of Mn and Fe solid and dotted lines show the corresponding trends observed in 
Milky Way clouds and their associated uncertainties (see text for discussion). 



consistent chemical evolution model is the subject of the 
following sections. 



3. DLAs as disk galaxies 

The details of our model for the chemical evolution of 
spiral disks are described in Boissier and Prantzos (2000, 
1999; hereafter BP2000 and BP1999, respectively), while 
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its application to the evolution of Zn in DLAs is made 
in Prantzos and Boissier (2000a, hereafter PB2000a). In 
Sect. 3.1 we briefly recall the main features of the model 
and its success in reproducing a large number of observa- 
tions concerning present day disks. In Sect. 3.2 we present 
in some detail the only two novel ingredients with respect 
to BP2000 and BP1999, namely the metallicity dependent 
yields of Woosley and Weaver (1995, hereafter WW1995) 
for intermediate mass elements, and the new yields of SNIa 
from of Iwamoto et al. (1999). In Sec. 3.3 we extend our 
study to the evolution of intermediate mass elements ob- 
served in DLAs. 

3.1. Description of the model 

A detailed model of the chemical and spectro-photometric 
evolution of the Milky Way is constructed in BP1999. 
The disk of our Galaxy is simulated as an ensemble of 
concentric rings gradually built up by infall of gas of 
primordial composition. Wc adopted a Star Formation 
Rate (SFR) suggested by the theory of density waves 
ip cx T} g ^ s V(R)/R (where V(R) is the circular velocity 
at radius R), which is in agreement with the observed 
current SFR profile in the Milky Way. The infall rate is 
exponentially declining, with a timescale t = 7 Gyr in the 
local disk as to reproduce the observed G-dwarf metal- 
licity distribution. Inner galactic zones are formed earlier 
than the outer disk (formation "inside-out"). The number 
of observables in our Galaxy, concerning the history of the 
solar neighborhood as well as the current profiles of var- 
ious quantities, is larger than the few "free" parameters 
involved in this model and testifies to its validity. 

After "calibrating" the model on the Milky Way (in 
fact, the SFR efficiency and the infall timescales as a 
function of local surface densities) we extended it to the 
study of local spirals (BP2000). We used simple "scaling 
laws", suggested by Mo, Mao and White (1998) in the 
framework of Cold Dark Matter scenarios of galaxy for- 
mation. Within this simplified approach, disks are char- 
acterised by two parameters: circular velocity Vc (de- 
termined by the mass of the dark halo), and spin pa- 
rameter A (measuring the specific angular momentum 
of the dark halo). For instance, the disk scalelcngth is 
Rd = Rd.MW * (Vc /Vc,mw) * (A/Amw), where the index 
MW refers to the corresponding value in the Milky Way 
(see BP2000 for details). 

Our grid of models covers the range 80 < Vc(km/s) < 
360 and 1/3 < \/\mw < 3, which corresponds to cur- 
rent High Surface Brightness disks. Our results reproduce 
quite satisfactorily a large number of observed proper- 
ties of present-day disks: sizes, surface brightness, Tully- 
Fisher relationship, colours, abundances, integrated spec- 
tra, etc... (BP2000). Also, the model reproduces the ob- 
served abundance gradients in disk galaxies (Prantzos and 
Boissier, 2000b) and the gas fraction, star formation effi- 
ciency and effective ages of a large sample of nearby spi- 



rals (Boissier et al. 2000). We notice that, according to our 
models, massive disks form quite rapidly (in a few Gyr) 
while it takes several Gyr for the small ones to be formed. 

In order to apply the model to the study of DLAs, the 
following statistical factors were taken into account: 

First, at redshift z = 0, a velocity distribution func- 
tion can be obtained by combining a Schcchtcr luminosity 
function and an observed Tully-Fisher relationship, as in 
Gonzalez et al. (2000): 



F v (V)dV = ^* 



exp 



dV 



(1) 



For the parameters of Equ. (1) we adopted the values of 
the fifth row of Table 4 in Gonzalez et al. (2000), corre- 
sponding to the velocity interval covered by our models; 
another choice would have small impact on the results of 
this paper. 

Results of numerical simulations (see e.g. Mo, Mao and 
White, 1998) give the distribution function of the spin 
parameter A: 



F\(X)dX = 



2iro-\ 



exp 



ln 2 (X/X) 



2ai 



dX 

T 



(2) 



with A = 0.05 and a\ = 0.5. 

Finally, the probability that a line of sight to a QSO 
intercepts a disc in the radius interval [R, R + dR] is 
proportional to the geometrical cross-section Fn(R)dR = 
2TtRdR/-KR 2 M (where Rm is the radius of the largest disc 
in our models), favouring the detection of the outer regions 
of the larger disks. 

At any time, or redshift, the average value of a quan- 
tity Q (for instance an abundance ratio) in our models is 
computed from: 



{<D) _ Sx Ivg Sr F*(F)Q(\ V c ,R) dR dV c dX 
{) JxfvJn F WdRdV c dX 



(3) 



where F = F(X,V C ,R) = F V (V C )F X (X)F R (R). We im- 
plicitly assume that those three functions are mutually 
independent and also time-invariable. This is of course an 
over-simplification, valid as long as mergers and intereac- 
tions do not modify too much Fy and F\ . This is probably 
a good approximation at low redshifts, but becomes more 
and more problematic at high redshifts (z > 3). 

The function $ in Equ. (3) plays the role of a "filter" . 
To obtain the true mean value over the whole ensemble 
of disks, one should take $(.F) = 1. However, Boisse et 
al. (1998) argued that observations of DLAs suffer from 
selection biases, as discussed in Sec. 2. Measurement of 
zinc abundances vs. HI column density show an anti- 
correlation, resulting from a bias against too large or too 
small values of T = [Zn/H]+log(N(HI)). 

To account for that bias, we apply Equ. (3) with ^(J 7 ) 
= 1 for 18.8 < T < 21 and elsewhere; this gives the 
average value of Q over those parts of our disks that have 
DLA-likc properties (at least, as far as Zn abundances and 
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HI column densities are concerned). As shown in PB2000a, 
this empirical "filter" can explain the weak evolution with 
rcdshift of [Zn/H], and the corresponding average value 
and observed scatter among DLAs (see also top-left pannel 
in Fig. 2 and discussion in Sec. 3.3). In the present work 
we extend the model to the study of other elements; the 
adopted yields are presented in the next section. 

3. 2. Yields of massive stars and SNIa 

An important ingredient in studies of galactic chemical 
evolution is the stellar yields of various elements. Mas- 
sive stars are the main producers of most heavy isotopes 
in the Universe. Elements up to Ca are mostly produced 
in such stars by hydrostatic burning, whereas iron peak 
elements are produced by the final supernova explosion 
(SNII), and by white dwarfs exploding in binary systems 
as SNIa. Most of He, C, N and minor CO isotopes, as well 
as s-nuclei come from intermediate mass stars(2— 8M Q ). 
Since we are mainly concerned with the abundances of a 
and iron-peak elements in DLAs, we consider no yields 
from intermediate mass stars in this work. 

We use the metallicity dependent yields of WW1995, 
which are given for stars of mass M = 12, 13, 15, 18, 
20, 22, 25, 30, 35 and 40 M and metallicities Z/Z = 0, 
10~ 4 , 10~ 2 , 10- 1 and 1. The WW1995 yields, folded with a 
stellar Initial Mass Function (IMF) lead to approximately 
solar abundance ratio of O/Fe (or a element), without the 
Fe contribution of SNIa. This lead Timmes et al. (1995b) 
to suggest that the Fe yields of WW1995 are probably 
overestimated. Following their suggestion, we adopted in 
our model half the nominal values for the WW1995 yields 
for iron peak elements (from Cr— Zn). Taking into account 
the uncertainties currently affecting those yields, such a 
reduction is not unreasonable. Our procedure allows to 
produce the observed evolution of O/Fc and a/Fe in the 
Milky Way halo, and does not alter the abundance ratios 
between iron peak elements (which are most commonly 
observed in DLA systems). 

To account for the additional source of Fe-peak ele- 
ments, required to explain the observed decline of O/Fe 
abundance ratio in the Milky Way disk (e.g. GP2000), 
we utilise the recent yields of SNIa from the exploding 
Chandrashekhar-mass CO white dwarf models W7 and 
W70 of Iwamoto et al. (1999). These are updated versions 
of the original W7 model of Thielemann et al. (1986), cal- 
culated for metallicities Z = Z Q (W7) and Z = (W70), 
respectively. These SNIa models lead to an oveproduction 
of Ni, which will also be clearly seen when we compare our 
model results of [Ni/Zn] with observations in low redshift 
DLAs. 

Finaly, notice that, due to uncertainties in the 
Zn yields of WW1995 (see Timmes et al. 1995b and 
PB2000a), we assume that Zn behaves always as Fe and 
we use then our model values of [Fe/H] as predictions for 
[Zn/H]. This is justified by the fact that observations in 



the Milky Way disk and halo show [Fe/Zn] = for a large 
range of metallicities. However, we shall question that as- 
sumption in Sec. 4.2. 

In recent papers we have shown how these metallic- 
ity dependent yields, combined with the model described 
in Sec. 3.1, lead to satisfactory results concerning the 
chemical history of the solar neighborhood (GP2000) and 
the current abundance gradients of intermediate mass el- 
ements in the the disk of the Galaxy (Hou et al. 2000). 
The use of metallicity dependent yields is a necessary con- 
dition for a consistent study of the abundance patterns in 
DLAs. 

3.3. Metallicity evolution vs observations 

We calculated the evolution of the abundances of all iso- 
topes between H and Zn in the different zones of our disk 
models. In this section we present results for eight ele- 
ments (Zn, Al, Si, S, Mn, Cr, Fe and Ni) that have been 
observed in DLAs. 

In Fig. 2 we present the absolute abundances of each 
element as a function of redshift. We display our results in 
two different ways. First, we show in long-dashed curves 
the extreme values, i.e. the highest and lowest metallicities 
at a given redshift, obviously corresponding to the inner- 
most zones of massive disks and the outermost zones of 
low mass disks, respectively (defined in PB2000a and here 
to be the zones at radial distance of 0.5 Rd and 5.5 Rd, re- 
spectively). The corresponding average value, taking into 
account various statistical factors (as discussed in Sec. 3.1) 
is given by Equ. 3 with ^(J 7 ) = 1 in all zones and is shown 
by a dotted curve. It can be seen that for all elements this 
"unfiltered" average value, as well as the extreme values, 
display strong evolution and are in disagreement with the 
data. 

The shaded area of Fig. 2 is obtained by adopting the 
"filter" suggested by Boisse et al. (1998) for Zn and al- 
ready introduced in PB2000a, i.e. by putting in Equ. (3) 
= 1 for 18.8 < T < 21, and = otherwise. 

The upper value of this filter (T <21) is maintained for 
all elements, since for larger values, obscuration of back- 
ground quasars inhibits DLA detection (at least accord- 
ing to our interpretation of the observed Zn/H vs Nhi 
anticorrelation, see PB2000a). However, the lower limit 
[T > 18) adopted for Zn corresponds to the minimal col- 
umn density of Zn atoms detectable in current surveys; 
it should not be the same for other elements, since it de- 
pends (among other) on which ionisation lines are used 
for abundance determination. We adopt then as a con- 
servative lower limit for all other elements the condition 
log(Njn) > 20.2, which is simply the definition of a DLA. 
In other terms, we consider only the zones of our models 
that have surface densities Sgos >2 M Q pc~ 2 . 

Applying this filter to our results (i.e. by excluding 
the very- metal rich and dense inner zones, as well as the 
low-density metal poor zones) we find a much narrower 
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range of abundances at a given redshift, as well as a milder 
abundance evolution than in the original models (the solid 
curve in Fig. 2 gives the average value in the "filtered" 
zones, i.e. those with log(Nni)> 20.2 and T < 21). In 
PB2000a we found that this "biased evolution" picture is 
quite compatible with observations of Zn in DLAs. Here 
we show that it is also compatible with all available data 
on DLA abundances: indeed, both the observed dispersion 
at a given redshift and the mild overall evolution are nicely 
reproduced. 

We find that in our models the average abundance 
varies at a rate d < [a/H] > jdz ~ —0.25 for a— elements 
and d < [Fe/H] > jdz ~ —0.33 for Fe in the redshift range 
z = 0.5—4. In the latter case, the larger rate is due to the 
introduction of a late source of Fe, namely SNIa. We note 
that in a recent work Savaglio e al. (1999), after correct- 
ing for depletion a large number of adundances in DLAs 
find that the evolution of the intrinsic DLA metallicity Z 
as a function of redshift can be described by dlogZ/dz=- 
0.30±0.06 for 0.5< z <4, i.e. a value not very different 
from ours. 

The scatter of absolute abundances at a given redshift 
in our models is attributed both to the different types 
of galactic disks considered here (massive ones evolving 
more rapidly than low mass ones) and to the differential 
evolution of the various regions within a given disk (inner 
zones evolving more rapidly than the outer ones). The 
obtained scatter, typically a factor of 30—50 at a given 
redshift, reproduces fairly well the observed one. We notice 
that in Fig. 2 we find a smaller scatter in the "filtered" 
values of Zn than in those of Fe, although we assumed 
that both elements evolve in the same way (i.e. [Zn/Fe] = 
0). This is due to the different lower filters used: T > 18.8 
for Zn and log(N(HI))>20.2 for Fe. 

The mild evolution and the large corresponding scat- 
ter of metal abundances in DLAs have been addressed in 
previous works in the field. Observational papers attribute 
the former to selection bias (e.g. Vladilo et al. 2000) al- 
though no quantitative model is proposed to substantiate 
this argument. On the other hand, most theoretical pa- 
pers (e.g. Timmes et al. 1996, Ferrini et al. 1997, Prantzos 
and Silk 1998, Jimenez et al. 1999, Lindner et al. 1999) ig- 
nore such biases and find invariably a metallicity evolution 
stronger than observed. Selection biases are properly con- 
sidered in PB2000a. Meusinger and Thon (1999) also take 
such effects into account, but consider only disk galax- 
ies in a narrow range of masses (5—8 10 10 M Q ), i.e. with 
quite similar star formation histories. As a result, even af- 
ter considering selection effects, they find that the average 
metallicity in their models increases by a factor of ~ 10 
between z = 1 and 3, much larger than suggested by the 
observations. For that same reason, they fail to reproduce 
the observed scatter at a given redshift and they obtain 
much smaller values. 
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Fig. 2. Evolution of elemental abundances in our disk 
galaxy models for 8 elements, Al, Si, S, Cr, Mn, Fe, Ni 
and Zn, as a function of redshift z. Observational data 
for DLAs are from references listed in Table 1, with typ- 
ical errors shown in the bottom left corner of each panel. 
In each panel, the shaded area includes all the "filtered" 
zones of our models (i.e. by taking various biases into ac- 
count, see text) and the thick curve is the corresponding 
mean value in those same zones. A weak evolution and an 
important scatter, both compatible with observations, are 
obtained for all elements. For comparison, we also give the 
results corresponding to the "unfiltered situation" (i.e. all 
the zones of our models): the long dashed curves corre- 
spond the upper and lower limits, while the short dashed 
curve is the corresponding mean value. Our disk models 
(which include no correction for dust depletion effects) 
have a true metallicity evolution stronger than the one 
observed in DLAs. 



4. Abundance ratios and depletion patterns in 
DLAs 

Abundance ratios are more reliable tracers of the chemical 
evolution than absolute abundances. They allow to iden- 
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tify whether an element originated from long lived sources 
(low mass stars or SNIa) or short lived ones (massive stars 
and SNII), or whether it is affected by the odd-even effect. 
The analysis of abundance ratios in the Milky Way disk 
and halo stars does provide important insight on the chem- 
ical history of the Milky Way as well as on the different 
production sites of various elements (despite the existence 
of large uncertainties in both theory and observations, see 
GP2000). 

However, situations in DLAs are more complicated. 
The abundances observed in DLAs may not represent 
the real chemical composition of the system if part of 
the elements is removed from the gas to the solid phase 
(dust grains), as happens in the interstellar medium of our 
Galaxy. In the discussion of Sec. 3.3, we did not consider 
any such kind of dust depletion effects. We have argued, 
in Sec. 2, that the observed abundance ratios [X/Zn] in 
DLAs, when plotted as a function of T, suggest that some 
depletion does occur in those systems, at least for Fe, Cr, 
Mn and Ni. We have also argued that the depletion pat- 
tern for Fe, when compared to the one observed in the 
Milky Way, suggests that the line of sight intercepts sev- 
eral gas clouds in a given DLA. 

In this section, we compare observed abundance pat- 
terns in DLAs with those obtained in our models of disk 
galaxy evolution. Our aim is to find out : i) whether there 
are any discrepancies between models and observations, 
ii) whether such discrepancies can be attributed to de- 
pletion, and iii) whether the resulting depletion patterns 
correspond to well known ones, such as those observed in 
the local disk or halo clouds. 

4-.1. Abundance ratios in DLAs 

In Fig. 3 we plot the evolution of abundance ratios [X/Zn] 
as a function of metallicity ( = [Zn/H]), redshift and metal 
column density ( T ) for the elements Al, Si, S, Cr, Mn, 
Fe and Ni. As in many other works, we assume that Zn is 
not depleted at all in dust and, therefore, it can be used as 
a tracer of metallicity. As explained in Sec. 3.2, because 
of uncertainties about the nucleosynthesis of Zn, we as- 
sume that [Fe/Zn] = always in our model. The heavy 
shaded regions in Fig. 3 present the results corresponding 
to the "filtered" zones of our models (i.e. satisfying DLA 
constraints), while the light shaded ones correspond to all 
other zones. 

a elements Si and S: a elements present a well known 
behaviour in our Milky Way disk and halo. The [a/Fe] 
ratio is nearly constant in the halo, at a value of 0.3 ~ 
0.5, and declines gradually in the disk from [Fe/H] <~ —1 
to 0. The later feature is attributed to the late (~ 1 Gyr) 
contribution of SNIa to the disk composition. 

The observed [Si/Zn] ratios in DLAs do not show any 
trend with metallicity, redshift or column density, and 
they are solar, on average (taking into account observa- 
tional uncertainties). Our model results, when plotted as 



a function of metallicity or column density cannot be di- 
rectly compared with the corresponding data, because a 
given metallicity in the model may be reached at widely 
different epochs, dependent on the chemical history of the 
corresponding region. Only when [X/Zn] is plotted as a 
function of redshift can a meaningful comparison to the 
data be made. From the middle panels of Fig. 3 it can be 
seen that at early times (z > 2.5) the model [Si/Zn] <~ 0.4, 
with a very small variation around that value. 

At late times [Si/Zn] declines in general, albeit with 
large variations. The decline is due to the late appearence 
of Fe producing SNIa and to our assumption that [Fe/Zn] 
= always. There are large variations of [Si/Zn] between 
the various zones of our models (light shaded area), be- 
cause SNIa enter the scene at different times depending 
on the previous SFR history of the corresponding zones. 
In the "filtered" zones, however, the scatter around the 
mean value is very small. Compared to that well defined 
model trend, the data show systematically lower values. 
This difference could, in principle, be attributed to deple- 
tion by ~ 0.4 dex, an amount comparable with the one 
seen in the local warm clouds of the disk and the halo. 

However, [S/Zn] in our models behaves exactly as 
[Si/Zn] (for the same reasons) and the corresponding DLA 
data suggest again the same degree of depletion. This is 
difficult to understand, since S is not expected to be de- 
pleted at all. Ockam's razor, along with the fact that S/Zn 
and Si/Zn in DLAs are independent of Zn column density 
(Fig. 1) suggest then that these two elements are not de- 
pleted in DLAs; they also suggest that the intrinsic DLA 
ratios of those elements are ~ solar, irrespectively of red- 
shift. 

This conclusion does not answer the question about 
the nature of DLAs. Several authors (Centurion ct al. 
2000; Jimenez et al. 1999, Vladilo 1999) suggested that 
Low Surface Brightness galaxies (LSBs) are good candi- 
dates, since they probably have different nuclcosynthctic 
histories (at present unknown and unconstrained). At this 
point we emphasize that our scenario involves many differ- 
ent chemical histories: disks of different masses form the 
bulk of their stars at different epochs and within them dif- 
ferent zones evolve at different rates. As a result, at late 
times, the [a/Zn] ratio varies considerably in our mod- 
els. But the zones satisfying the DLA criteria ( T > 18.8 
for Zn and log(N(HI))> 20.2 for other elements) display 
a narrow range of values at all times. We shall return to 
that point in Sec. 4.2. 

Fe peak Elements: Cr, Mn, Fe and Ni are produced 
by massive stars in the early life of the Milky Way and 
mostly by SNIa at late times. Observations of their abun- 
dances in halo and disk stars show that Cr, Fe and Ni 
behave similarly, at least down to [Fe/H]~ —2, while the 
monoisotopic Mn behaves differently: the [Mn/Fe] ratio is 
~-0.4 in the halo and increases to its solar value in the 
disk. This behaviour is well accounted for by the metal- 
licity dependent yields of WW1995 for massive stars and 
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Fig. 3. Abundance ratios versus metallicity (left) redshift (middle) and metal column density T (right). In all panels, 
the light shaded area includes all the zones of our models, while the heavy shaded area only those corresponding to 
the "filtered" zones (i.e. satisfying DLA constraints, see Sec. 3.3). Note that [Fe/Zn] = in our models, because we 
have assumed that [Zn/H] = [Fe/H] always (see Sec. 3.2). Data for DLAs are from Table 1 and typical error bars 
are displayed in the middle panels. Deviation of the data from the model values at a given redshift is, in principle, 
attributed to dust depletion, which is not included in the models (see discussion in Sec. 4.2). 



those of Iwamoto et al. (1999) for SNIa, as shown by de- 
tailed chemical evolution models in GP2000. 

Our models show that the [Cr/Zn] and (by assump- 
tion) [Fe/Zn] ratios are ~ solar and independent of metal- 
licity or redshift. As expected, [Mn/Fe]~ —0.4 at low 



metallicities (high redshifts) and increases slowly at late 
epochs. The small overproduction obtained for [Ni/Zn] ob- 
tained at late times has a well known origin: the adopted 
SNIa yields of Iwamoto et al. (1999) which oveproduce Ni 
by a factor of ~ 2 (see GP2000 for a discussion). Taking 
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[Ni/Zn] = over the whole metallicity range would thus 
be more appropriate. 

In all four cases, the observed [X/Zn] ratio in DLAs is 
clearly below solar or model values, pointing to depletion 
as the most obvious explanation. This conclusion is cor- 
roborated by the observed anticorrelation of [X/Zn] with 
Zn column density, as argued in Sec. 2. 

In a recent work Pettini et al. (2000) selected a small 
subsample of [Mn/Fe] values in DLAs with [Zn/Cr] < 0.3, 
i.e. with small amount of depletion; assuming then that 
all Fe peak elements are depleted by the same (small) 
amount, they found that [Mn/Fe]~ —0.4 in DLAs, inde- 
pendent of metallicity. This picture is not quite compatible 
with the one in the Milky Way, but because of their small 
sample it is difficult to conclude whether their finding has 
implications for the role of SNIa in DLAs or not. 

4-2. Depletion patterns 

From the discussion of the previous section it turns out 
that several elements in DLAs are depleted, at least with 
respect to the predictions of our chemical evolution mod- 
els, which reproduce satisfactory several DLA features 
(e.g. Fig. 2). It is interesting to see how the predicted 
depletion patterns compare to some well known ones, as 
those of the warm halo and disk clouds in the Milky Way. 

In the upper panels of Fig. 4, we compare the observed 
[X/Zn] ratios in DLAs to the corresponding ones of our 
models i.e. at the absorbers' redshifts, (as given in the mid- 
dle panels of Fig. 3). Taking into account the discussion 
of Sec. 2, namely the fact that the depletion patterns de- 
pend on Zn column density, we split the data into low 
and high T values [T < 19.8 and T > 19.8, respectively). 
We also show, for comparison (in both panels) the deple- 
tion patterns in the Milky Way warm halo and disk clouds 
(open symbols, connected by solid lines, data from Savage 
& Sembach 1996). Due to the uncertainties (unknown in- 
trinsic abundances in DLAs, measurement errors), such a 
comparison has no sense for individual DLA, but only in 
a statistical sense. 

In the left upper panel of Fig. 4, it is seen that the de- 
pletion amount obtained at low T is, on average, smaller 
than in warm halo clouds. The average depletion is sim- 
ilar for all elements, <~ —0.2 dex, and comparable to the 
error bars of the measurements. In the right upper panel 
of Fig. 4, the depiction is slightly larger for all elements. 
The depletion pattern ressembles, on average, the one of 
warm halo clouds for Fe, Ni and Cr, but also for Mn and 
Si (i.e. within error bars). However, as dicussed in the pre- 
vious section, it is simply impossible to accomodate S in 
this scheme, since it is hard to imagine that its ratio to 
Zn can be subsolar. 

The problem with S may be due to the fact that we use 
as reference pattern the one resulting from our models. We 
checked then whether using the solar abundance ratios as 
reference pattern, could lead to a better fit of the data. 



The results of that comparison appear in the lower panels 
of Fig. 4, again for low and high T values (left and right 
lower panels, respectively). At low T hardly any depletion 
at all is found, for all elements except Mn; this is puzzling, 
since Mn can be no more depleted than the other Fe peak 
elements. At high T, the resulting depletion pattern com- 
pares fairly well with the one of the warm halo clouds for 
all elements. 

In summary: if one assumes that the intrinsic metal- 
licity pattern in DLAs is the one given by our models, a 
difficulty arises with S which is found to be depleted at 
high T . If on the other hand it is assumed that the intrin- 
sic DLA pattern is solar, a problem arises with Mn, which 
is found to be depleted at low T. 

We note that our assumption of [Fe/Zn]=0 always, al- 
though based on observations of Milky Way stars, may be 
inconsistent and potentially wrong. Indeed, it implies that 
^2/3 of Zn are produced at late times by SNIa, exactly 
as is Fe in the solar vicinity (see GP2000). This assump- 
tion introduces a time-delay in the production of the bulk 
of Zn with respect to the one of other a-elements, like S, 
coming from massive stars. The reason for that delay is 
that the bulk of Fe is always produced with some delay 
in our models, through SNIa explosions. This explains the 
obtained decline of a/Zn in our models only at redshifts 
lower than However, physically, this picture is not cor- 
rect: it is well known that SNIa do not produce Zn (e.g. 
Iwamoto et al. 1999). The fact that in the solar vicinity 
Zn behaves as Fe also at late times is obviously a coinci- 
dence. The late Fe production of SNIa is matched by the 
Zn production of SNII, which is metallicity dependent in 
both the yields of WW1995 as well as in the recent work 
of Limongi et al. (2000). The yields of WW1995, however, 
do not reproduce the observed evolution of Zn/Fe in the 
Milky Way and we used Fe as a proxy for Zn; in that way 
we "forced" a pure metallicity effect (which can, in princi- 
ple, manifest itself quite early in zones that evolve rapidly 
to high metallicitics) to mimic a time-delay effect. In that 
sense, our treatment of Zn is not quite correct. Obviously, 
the nucleosynthesis of Zn should be carefully scrutinised 
in future models of massive star nucleosynthesis. 

At this point it seems difficult to gain more insight into 
the nature of DLAs by simply considering their abundance 
patterns. We think, however, that progress can be made 
through further observations of S in DLAs of high Zn col- 
umn density and of Mn in DLAs of low Zn column den- 
sity; such observations would certainly help to establish 
whether the aforementioned trends are real or not. 

5. Summary 

In this work we study the chemical abundance patterns 
observed in DLAs. First, we show that there is strong 
observational evidence for depletion of the Fe— peak ele- 
ments, obtained when their abundance ratios [X/Zn] are 
plotted as a function of Zn column density. The resulting 
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Fig. 4. Depletion patterns in DLAs. In the upper panels, observed DLA abundance ratios [X/Zn] are compared to 
those resulting from our models, taken at the absorber's redshift. In the lower panels, the intrinsic DLA abundance 
pattern is simply assumed to be solar. In both cases, the DLA data are plotted according to the corresponding Zn 
column density, given by T = [Zn/H]+log(NHi): left panels correspond to T < 19.8 and right panels to T > 19.8. In 
all panels the depletion patterns obtained in warm halo and disk clouds in the Milky Way are shown by open symbols, 
connected with solid lines (data from Savage and Sembach 1996). Typical error bars for measured [X/Zn] ratios in 
DLAs are displayed in the right upper panel. 



anticorrelation is reminiscent of similar depletion patterns 
obtained in Milky Way clouds as a function of gas col- 
umn density. However, the intrinsic gas metallicity does 
not differ much from solar in Milky Way clouds, while it 
varies by factors ~ 100 in DLAs. For that reason we argue 
that Zn column density is a more appropriate parameter 
to measure depletion, since it is proportional (to a first 
approximation) to the amount of dust along the line of 
sight. Further measurements of [X/Zn] ratios as a func- 
tion of T = [Zn/H]+log(N H i) in DLAs will help to put 
the anticorrelation found here on a firmer basis. 

Assuming that DLAs are (proto-)galactic disks, we 
study their chemical evolution with detailed models, mak- 
ing use of metallicity-dependent yields from massive stars 
and SNIa; these models reproduce quite successfully a 
large number of properties of spirals at low and high 
redshift as shown in a series of papers (Boissier et al. 
2000; BP2000, 1999; PB2000a, 2000b). By taking into ac- 
count observational constraints on the gas and metal col- 
umn densities of DLAs (namely that log(Nni> 20.2 and 
18. 8< T < 21) we find that our models reproduce fairly 



well (and without including any depletion effects) the ob- 
served weak evolution of DLA metallicity with redshift 
and the corresponding scatter, for several elements de- 
tected in DLAs: S, Si, Mn, Cr, Fe, Ni and Zn. 

Encouraged by the success of the model we study the 
corresponding predicted abundance patterns. We find no 
clear correlation of [X/Zn] with metallicity or HI column 
density. We also find that Si/Zn and S/Zn are supersolar 
at high redshifts, while Mn/Zn is subsolar. The latter ef- 
fect results from the metallicity dependent yields adopted 
here and is also observed in the Milky Way (see GP2000 
and references therein). The former results from our as- 
sumption that [Zn/Fe] = always, assumption based on 
observations in the Milky Way and made because of the 
poorly understanding of Zn nucleosynthesis. 

Comparing the observations of [X/Zn] in DLAs to our 
model calculations ((which do not include the effects of 
dust depletion) we find evidence for dust depletion effects 
in the observations for Cr, Ni, Fe, Mn and Si, for high T > 
19.8. This depletion is ~0.4 dex, on average, and compa- 
rable to the one observed in warm halo clouds. However, 
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a difficulty arises with S which is found to be depleted at 
high T . If on the other hand it is assumed that the intrin- 
sic DLA pattern is solar, a problem arises with Mn, which 
is found to be depleted at low T. 

We think that currently available data on DLA abun- 
dances do not allow to establish the intrinsic abundance 
patterns of those systems. Progress can be made through 
further observations of abundances in DLAs, in particular 
those of Mn and S, by taking into account the aforemen- 
tionned role of the metal column density. We also note 
that the conclusions of this work are largely based on the, 
presently poorly understood, nucleosynthesis of Zn; the 
production of that element should be carefully studied in 
new models of massive star nucleosynthesis. 
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